,  i  i 


AD-A278  459 

— — 

OFFICE  OF  NAVAL  RESEARCH 
Grant  N00014-93-1-0351 
RAT  Code  413302Q — 09 

Technical  Report  No.7 

The  Kinetics  of  Heterogeneous  Electron  Transfer  Reactions  in  Polar  Solvents 

by 

W.  Ronald  Fawcett  and  Martin  Opallo 

Prepared  for  Publication 
in 

Angewandte  Chemic  International  Edition  in  English 

Department  Of  Chemistry 
University  of  California 
Davis,  CA  95616 


April  20, 1994 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government 


"This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited" 


94-12409 


euc  rr?  tctt'scted 3 


94  4 


o*, 

Au'  Am’ 


1 


A 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No  0704-0188 

-i.-ti# - *  3J»d  »-3  •♦OitiF'J  »««g  .  '•'JB'g.rq  !f  %  Ol*'a**  O'  in.  stotet  o*  ?»!.* 

*: ■  ,  .  .  ‘ _ _  -  •  .3^;w4-aest  ->  •  *  *3wv  7*3  •-■ s  ourog*  ':  ••tNg^irtrs  Oo*Mt'on$  sno  *#ocm  f2*»  e^erso" 

,f  s  <  t-j  ,  *;•;-  i  ji.*;  j^-g  •;  '•  3*  Ml"W*’"T  ♦'*0  3  -0^0*  a  »0«',n»CF<  O'*  *"!iect  3  K*  -0  ’  38)  A**P  nqtO"  -C*350J 

T‘  AGENCY  USE  ONLY  ([five  Qitnk)  2  REPORT  OATE  3.  REPORT  TYPE  ANO  OATES  COVEREO 

April  20,  1994  Technical 

4.  TITLE  ANO  SUBTITLE 

The  Kinetics  of  Heterogeneous  Electron  Transfer 

Reactions  in  Polar  Solvents. 

5.  FUNDING  NUMBERS 

N00014-93-1-0351 

6.  AUTHOR(S) 

W.R.  Fawcett  and  M.  Opallo 

7.  PERFORMING  ORGANIZATION  NAM£(S)  ANO  AOORESS(ES) 

Department  of  Chemistry 

University  of  California 

Davis,  CA  95616 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

No.  7 

9.  SPONSORING' MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy 

Arlington,  VA  22217-5000 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

Prepared  for  publication  in  Angewandte  Chemie  International  Edition  in 

English. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Unclassified 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Data  for  heterogeneous  electron  transfer  reactions  obtained  in  the  last  ten  years  have  been  examined 
within  the  context  of  contemporary  theory.  Special  attention  is  focussed  on  systems  for  which  rate 
constants  and  activation  parameters  are  available  as  a  function  of  the  solvent,  and  as  a  function  of 
temperature.  The  role  of  reactant  structure  in  determining  the  kinetic  parameters  is  also  considered. 

Double  layer  effects  both  at  unmodified  and  modified  electrodes  are  discussed  in  detail. 

Experimental  techniques,  especially  those  suited  to  measuring  fast  reactions,  are  also  outlined. 

14.  SUBJECT  TERMS 

15.  NUR8BER  OP  PAGES 

IB.  PRICE  COOE 

17.  SECURITY  CLASSIFICATION 

Of  REPORT 
unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

unclassified 

19.  SECURITY  CLASSIFICATION 
OP  ABSTRACT 

unclassified 

20.  UMTATION  OP  ABSTRACT 

unclassified 

V>N  7540-01  -280-5500  Standard  Form  191  (Rev  2-89) 


nrnnwi  a*  amt  sta  *».n 
m-<s2 


1 


A  Review 


The  Kinetics  of  Heterogeneous  Electron  Transfer  Reactions  in 

Polar  Solvents 

by  William  Ronald  Fawcett41  and  Marcin  Opallo 


Data  for  heterogeneous  electron  transfer  reactions  obtained  in  die  last  ten  yean  have 
been  examined  within  the  context  of  contemporary  theory.  Special  attention  is  focused  on 
systems  for  which  rate  constants  and  activation  parameters  are  available  as  a  function  of  the 
solvent,  and  as  a  function  of  temperature.  The  role  of  reactant  structure  in  determining  the 
kinetic  parameters  is  also  considered.  Double  layer  effects  both  at  unmodified  and 
modified  electrodes  are  discussed  in  detail.  Experimental  techniques,  especially  those 
suited  to  measuring  fast  reactions,  are  also  outlined. 
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I .  Introduction 

Since  the  early  1980's  there  has  been  a  great  increase  in  interest  in  medium  effects  for 
election  transfer  (ET)  reactions.  This  revival  was  spurred  by  the  appearance  of  several 
theoretical  papers  [1-4]  which  predicted  the  existence  of  solvent  dynamical  effects  in  the 
kinetics  of  ET.  The  first  experimental  evidence  for  a  dynamical  efFect  in  a  heterogeneous  ET 
reaction  was  presented  by  Kapturkiewicz  and  Be  hr  [3].  This  was  quickly  followed  by  a 
report  describing  similar  effects  for  a  homogeneous  ET  reaction  [6].  Since  then  a  significant 
body  of  experimental  evidence  has  been  described  in  the  literature  relating  to  solvent  and 
other  medium  effects.  These  experiments  have  focused  on  three  areas  namely,  homogeneous 
self-exchange  ET,  heterogeneous  ET,  and  time-resolved  fluorescence  spectroscopy  involving 
intramolecular  ET.  At  the  same  time,  significant  advances  have  been  made  in  die  theoretical 
description  of  ET  processes. 

Several  recent  reviews  have  focussed  on  the  development  of  theory  of  this  elementary 
reaction  and  its  application  to  experimental  results  [7-9].  Heitele  [9]  has  examined  in  a 
particularly  lucid  and  useful  manner  the  theory  of  dynamic  solvent  effects  in  ET.  On  the 
experimental  side,  Barbara  and  Jarzeba  [10]  have  reviewed  experimental  results  in  time 
resolved  fluorescence  spectroscopy.  Weaver  [1 1]  has  discussed  data  for  homogeneous  ET 
reactions,  especially  those  obtained  in  his  laboratory. 

The  subject  of  the  present  review  is  heterogeneous  ET.  The  discussion  is  limited  to 
simple  reactions,  that  is,  those  in  which  both  oxidized  and  reduced  forms  remain  in  solution, 
only  one  election  is  transferred,  and  no  chemical  bonds  are  broken  or  formed.  Since  the 
original  work  of  Kapturkiewicz  and  Behr  [S]  on  dynamic  solvent  effects,  a  large  number  of 
related  studies  have  appeared  in  the  literature.  Our  goal  is  to  examine  these  results  within  the 
context  of  contemporary  theoretical  development  and  show  their  relationship  to  other 
experiments  describing  medium  effects  in  ET  No  attempt  is  made  here  to  discuss  *heory  in 
detail  since  this  has  already  been  done  [7-9].  However,  a  few  aspects  of  theory  which  were 


not  considered  in  previous  reviews  and  which  are  particularly  relevant  to  the  analysis  of 
experimental  data  are  considered.  As  well  as  presenting  the  important  results  obtained, 
advances  in  experimental  techniques  are  also  discussed.  It  is  important  to  emphasize  that  the 
improved  understanding  of  the  role  of  the  medium  in  ET  has  been  due  in  great  part  to  the  fact 
that  faster  processes  in  condensed  media  can  now  be  followed  experimentally  with  precision. 
Finally,  it  should  be  pointed  out  that  heterogeneous  ET  is  usually  studied  in  the  presence  of 
an  inert  electrolyte  in  order  to  improve  solution  conductivity  and  maintain  diffusion  as  a 
means  of  mass  transfer.  As  a  result,  the  medium  for  the  ET  reaction  involves  not  only  the 
solvent  but  the  ions  of  the  electrolyte  which  form  the  double  layer  at  the  electrode.  Double 
layer  effects  play  an  important  role  in  heterogeneous  ET  and  are  also  considered  in  this 
review. 

Theory 

2 . 1  Estimation  of  the  Activation  Parameters 

One  of  the  important  problems  considered  in  recent  years  is  improvement  of  the 

estimation  of  the  Gibbs  energy  of  activation  AG*  and  related  kinetic  parameters.  AO*  is 

* 

traditionally  regarded  as  composed  of  two  parts,  the  inner  sphere  contribution  AG^  and  the 

* 

outer  sphere  contribution  AGQS  [  12].  Estimation  of  the  former  requires  a  description  of  the 

changes  in  bond  lengths  and  angles  which  accompany  the  activation  process,  and  is  often 
based  on  the  appropriate  spectroscopic  data  [13].  The  outer  sphere  component  which 
accounts  for  the  work  done  to  reorganize  the  solvent  atmosphere  around  the  reactant  has 
usually  been  estimated  on  the  basis  of  a  model  based  on  the  Bom  theory  of  ionic  solvation 
[12, 13].  In  the  case  of  the  simple  heterogeneous  ET  reaction, 

A  +  e  — >  B  (1) 


this  contribution  is  given  by 
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...  Nq5q_  n  J_\  /J_  In 

32hEq  lr*  RJ  (eop  ej 


(2) 


where  &>  is  the  Avogadro  constant.  So,  the  fundamental  electronic  charge,  Ia»  the  radius  of 

the  reactant  A  represented  as  a  sphere,  R>  the  distance  from  its  image  in  the  electrode,  £op,  the 

opdcal  relative  permittivity  of  the  solvent,  £$,  the  static  value,  and  £0,  the  permittivity  of  free 

space.  An  important  point  in  application  of  this  relationship  is  estimation  of  Ia  since  most 

reactants  are  not  spherical.  This  aspect  is  considered  in  mem  detail  below.  In  estimating 
* 

AGos  in  recent  work,  it  has  usually  been  assumed  that  die  distance,  R,  is  effectively  infinity, 
the  stabilizing  effect  of  any  image  formed  in  the  electrode  being  neglected.  Certainly,  in 
cases  where  the  ET  reaction  remains  simple  in  the  sense  that  the  reactant  and  product  do  not 
interact  with  the  electrode  itself,  R  is  equal  to  2i*  +  2di  where  &  is  the  thickness  of  die  layer 
composed  of  solvent  and  closest  approaching  counter  ions  between  the  reactant  and  die 
geometrical  electrode/solution  interface.  Thus,  R  corresponds  at  least  to  three  molecular 
diameters,  that  is,  it  is  the  order  of  six  times  Ia-  Moreover,  the  efficiency  with  which  the 
image  is  formed  in  the  electrode  can  be  questioned  on  the  basis  of  the  spatial  dispersion  of  the 

solvent  dielectric  permittivity  and  electric  field  penetration  into  the  metal  electrode  [14,  IS]. 

* 

Komyshev  et  al.  [14,  IS]  applied  the  non-local  electrostatic  treatment  to  estimation  of  AGQS 
for  both  heterogeneous  and  homogeneous  ET  reactions,  and  concluded  that  this  quantity 
actually  increases  as  the  reaction  site  moves  closer  to  the  electrode  instead  of  decreasing  as 
predicted  by  equation  (2).  However,  the  sharp  boundary  model  used  in  these  calculations 
[  IS]  can  be  criticized  because  it  is  not  applicable  in  regions  where  the  dielectric  properties  of 
the  system  are  changing  markedly  over  short  distances,  that  is,  near  a  reactant  at  the  reaction 
site  in  the  double  layer.  Nevertheless,  on  the  basis  of  a  comparison  of  activation  enthalpies 

for  the  same  reaction  under  homogeneous  and  heterogeneous  conditions  (vide  infra)  it  s.  ;ms 

* 

clear  that  R  is  sufficiently  large  that  1/R  can  be  neglected  in  estimating  AGQ5  for  the 
heterogeneous  process.  This  assumption  is  made  in  the  following  discussion. 
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Further  improvements  in  the  estimation  of  AGQ.  have  involved  the  use  of  non- 

primitive  models  for  ionic  solvation,  especially  the  mean  spherical  approximation  (MSA)  [16, 

17].  On  the  basis  of  the  hypothetical  two  step  charging  process  associated  with  thermal 

* 

electron  transfer  described  by  Marcus  [18],  AGQ.  can  be  regarded  as  giving  the  difference 
between  the  Gibbs  solvation  energies  of  the  reacting  system  at  optical  frequencies  and  at 
static  frequency.  This  can  be  seen  when  equation  (2)  is  rewritten  as 


327C£oTa 


/,  1  N  ,  NoCo  /. 

.  1\ 

(  Cop)  32nEorA  ( 

Cs) 

(3) 


in  a  form  which  clearly  shows  the  relationship  of  the  Marcus  model  to  the  Bom  continuum 

model  for  ionic  solvation.  Following  earlier  work  of  Wolynes  [19],  McManis  and  Weaver 

* 

[20]  wrote  expression  (4)  for  AGps  in  the  MSA: 


AG^s  =  -  [Yl  -  — )  - —— 

32jcEo  Lv  Copy  (rA+rsAop) 


(fA+rs/Xs). 


(4) 


where  is  is  the  radius  of  a  solvent  molecule  represented  as  a  sphere,  and  Aop  and  As, 
polarization  parameters  estimated  from  £op  and  £s,  respectively  by  equations  (5)  and  (6). 


16  Cop  -  Xgp(l  +  Xop)4 

(5) 

16  Es  =  Xjd+Xs)4 

(6) 

In  its  application  to  ion  solvation  at  static  frequencies  [21],  die  parameter.  A* 

accounts  for  the  polarizability,  softness,  and  non-sphericity  of  the  solvent  molecule  which  is 

otherwise  represented  as  a  hard  sphere  with  a  central  point  dipole.  For  polar  solvents,  the 

value  of  As  is  in  the  range  2.0  to  2.9  [22],  so  that  the  MSA  estimate  of  the  static  contribution 
* 

'oAQos  is  considerably  reduced  from  the  Bom  estimate.  In  the  case  of  the  high  frequency 
solvation,  Acp  is  a  number  close  to  unity,  so  that  the  estimate  of  the  optical  contribution  to 
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* 

AGos  is  also  greatly  reduced  in  magnitude  from  the  continuum  value.  However,  there  is  a 
serious  question  whether  the  MSA  model  developed  for  a  system  of  hard  spheres  with  point 
charges  and  point  dipoles  at  static  frequencies  [19-22]  can  be  applied  at  optical  frequencies 
where  there  is  no  c  ientational  contribution  to  dipolar  polarization  [23],  Thus,  the  McManis- 
Weaver  (MW)  model  [20]  can  be  questioned  within  the  context  of  the  MSA. 

* 

A  quite  different  approach  to  applying  the  MSA  to  estimation  of  AGQS  was  taken  by 

Fawcett  and  Blum  (FB)  [23].  Firstly,  they  assumed  that  the  optical  contribution  is  given  to  a 

good  approximation  by  the  Bom  estimate.  Secondly,  the  polarization  parameter  at  static 

frequencies  was  obtained  from  estimates  of  the  Gibbs  solvation  energies  of  simple 

monoatomic  monovalent  ions  rather  than  from  equation  (6).  As  a  result,  the  expression  for 
« 

AGQ5  becomes 


NqCj  f/ 

1  .  JL\  1  . 

<X  l\  1  1 

32xeo  U 

£op)  fA  1 

<  ej  (rA+rsAi)J 

(7) 


where  the  parameter  Xi  depends  on  whether  the  ion  is  a  cation  or  an  anion  but  not  on  its  size. 
Stated  otherwise,  the  FB  model  predicts  that  the  Gibbs  activation  energy  to  form  a 
monovalent  cation  from  a  neutral  molecule  via  electrooxidation  is  different  from  that  to  form  a 
monovalent  anion  from  the  same  reactant  via  electroreduction,  all  other  system  parameters 
remaining  the  same.  This  follows  from  the  experimental  observation  that  the  extent  of 
solvation  of  cations  and  anions  of  the  same  radius  in  a  given  solvent  is  not  the  same  [24]. 

* 

Two  other  activation  parameters,  namely,  the  enthalpy  AHQ£  and  entropy  of 
* 

activation  AS^  may  be  estimated  given  the  temperature  coefficient  of  die  Gibbs  energy. 

* 

From  equation  (3)  estimation  of  ASps  on  the  basis  of  the  Bom  model  requires  knowledge  of 
the  temperature  coefficients  d£g/dT  and  d£np/dT  which  are  available  for  all  polar  solvents  used 
in  electrochemical  studies  [24, 25].  In  the  case  of  the  MW  model,  one  also  needs  the 
temperature  coefficients  dX^p/dT  and  dX«,/dT  which  may  be  obtained  from  the  experimental 
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coefficients  d£«/dT  and  de^p/dT  after  differentiating  equations  (5)  and  (6)  with  respect  to 

temperature  [22].  Finally,  one  needs  the  temperature  coefficient  dXi/dT  in  the  FB  model. 

These  quantities  were  extracted  from  experimental  solvation  entropy  data  for  the  alkali  metal 

*  * 
cations  and  halide  anions  [22].  Once  one  has  estimated  ASqS  » the  activation  enthalpy  AHqj 

*  * 

can  be  found  by  adding  to  the  estimate  of  AGQ£. 

*  41 

Values  of  AGQ£  and  AHQ£  estimated  by  the  various  models  described  here  in  twelve 

polar  solvents  are  summarized  in  Table  1 .  The  reactant  size  parameter  £a  was  set  equal  to 

0.37  nm,  a  value  appropriate  for  the  metallocenes  [26]  whose  ET  kinetics  have  been  studied 

extensively  both  homogeneously  and  heterogeneously  [11].  The  first  observation  one  may 

*  * 

make  is  that  the  estimates  of  AGQ£  and  AHQ£  differ  by  only  a  few  kJ  mol*1  at  the  most,  the 

difference  being  negligible  on  the  basis  of  the  MSA  for  most  solvents.  However,  estimaaes 

based  on  the  MSA  differ  considerably  depending  on  the  assumption  made  regarding  the  high 

frequency  contribution,  and  the  method  of  estimating  the  polarization  parameter  at  low 

* 

frequencies.  Estimates  of  AGQ£  on  the  basis  of  the  MW  model  are  about  90  percent  of  those 

obtained  by  the  continuum  Marcus  model.  On  the  other  hand,  using  the  FB  model,  die 
* 

estimates  of  AGQ£  fall  to  62  percent  of  the  Marcus  estimate  for  a  reaction  involving  cations, 
and  80  percent  for  a  reaction  involving  anions.  These  differences  are  significant  and  should 
be  seen  experimentally. 

The  major  problem  with  the  MSA  approach  is  that  a  proper  estimate  of  ion-dipole 
interactions  at  high  frequencies  has  yet  to  be  made.  Some  support  for  the  FB  assumption  is 
found  from  the  finding  that  the  MSA  estimates  of  the  dielectric  properties  of  a  polar  solvent  at 
optical  frequencies  agree  well  with  those  based  on  classical  models  [27].  However,  the 
separate  problem  of  ion-dipole  interactions  at  high  frequencies  needs  to  be  addressed  as  well. 
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2.2  Estimation  of  the  Pre-Exponential  Factor 

Much  has  been  written  already  about  the  pre-exponential  factor  in  die  rate  constant  far 
ET  especially  with  respect  to  dynamic  solvent  effects  [1-4, 7, 9, 1 1].  The  discussion  here  is 
limited  to  those  aspects  which  relate  to  heterogeneous  ET  and  to  appropriate  expressions  for 
the  rate  constant.  On  the  basis  of  transition  state  theory  formalism  [13],  the  standard 
heterogeneous  rate  constant  is  usually  written  as 

ks  =  kKp  vn  exp(-AG*/RT)  (8) 


where  K  is  the  electronic  transmission  coefficient,  j£p,  the  encounter  pre-equilibrium 
constant,  and  the  nuclear  barrier-crossing  frequency.  &  is  unity  for  an  adiabatic  reaction 
and  falls  below  unity  as  the  extent  of  diabaticity  increases.  The  extent  of  diabaticity  is 
connected  with  the  degree  of  interaction  between  the  free  energy  surfaces  for  the  reactant  and 
product  and  is  determined  by  the  electronic  matrix  coupling  element  Hw,  [13].  The  value  of 
Hah  can  also  affect  the  nuclear  barrier  crossing  frequency,  ¥n»  so  that  when  Hab  is  large  and 
the  barrier  rounded,  ¥n  falls  below  the  value  which  it  has  for  small  values  of  H*  [4, 28]. 
Simple  relationships  connecting  these  quantities  are  given  by  the  Landau-Zener  model  [13] 


according  to  which  the  electron  tunnelling  frequency  is  given  by  equation  (9), 

2H*  (  4n3  V/2 


vel  = 


No*1  l^AG’RT  J 


(9) 


where  Mo  is  Avogadro's  constant  and  h,  Planck's  constant  The  electronic  transmission 
coefficient  is  then  given  by  equation  (10). 


2[l-exp(-vei/2vn)] 

2-exp(-vei/2vn) 


GO) 


When  £  is  unity  and  the  reaction  is  adiabatic.  However,  when  &1  «  ¥n»  the 

expression  for  £  becomes 


K  =  vei/vn 


(ID 
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and  electron  tunnelling  controls  the  magnitude  of  the  pre-exponential  factor.  This  diabatic 
limit  occurs  in  cases  where  the  extent  of  interaction  between  the  reactant  and  electrode  is 
minimized,  for  instance,  in  the  presence  of  a  blocking  film  on  the  electrode.  In  this  case  H*b 
is  small,  ¥ei  falls,  and  *  is  given  by  equation  (11).  Obviously,  the  structure  of  the 
electrode/solution  interface  plays  an  important  role  in  determining  these  quantities. 


The  encounter  pre-equilibrium  constant  £p  has  been  defined  within  the  context  of  the 
encounter  pre-equilibrium  model  for  heterogeneous  ET  reactions  [29]  and  contains  any 
double  layer  terms  and  steric  effects.  In  the  absence  of  the  latter  effects,  it  has  been  estimated 
to  be  approximately  60  pm  for  a  spherical  reactant  [29].  Thus,  after  correction  for  double 
layer  effects,  the  upper  limit  tor  the  product  xKp  is  ~60  pm.  Lower  values  can  be  due  to 
steric  effects  or  departure  from  reaction  adiabaticity.  Because  of  the  very  approximate  nature 
of  the  models  used  to  estimate  ]£  and  j£p  individually,  it  is  preferable  to  cite  only  die  product 
icKp  on  the  basis  of  an  analysis  of  experimental  data  [30]. 


The  nuclear  barrier-crossing  frequency  depends  on  vibrational  modes  within  the 
reactant  and  solvent  modes  to  an  extent  which  is  contingent  upon  the  relative  importance  of 
the  inner  and  outer  sphere  contributions  to  the  activation  barrier.  When  the  reaction  is 
adiabatic  and  the  inner  sphere  contribution  to  AG*  is  small,  this  frequency  is  given  by 
equation  (12)  [1-4] 


vn  = 


tl/1 


'  AG*  1/2 
4tc  RT J 


(12) 


where  j|„  is  the  longitudinal  solvent  relaxation  time.  For  Debye  solvents  in  the  continuum 
approximation,  Zl  is  given  by  [31] 


tL 


EooTp 

£s 


(13) 


where  is  the  high  frequency  relative  solvent  permittivity  and  Id,  the  Debye  relaxation  time. 
Values  of  II  have  been  collected  for  a  large  number  of  polar  solvents  [32]  whose  behaviour 
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follows  the  Debye  model  to  a  good  approximation.  Thus,  the  pre -exponential  factor  is 
predicted  to  depend  on  the  solvent's  dielectric  relaxation  properties,  parameters  which  are 
quite  different  in  origin  to  those  which  determine  the  Gibbs  activation  energy  barrier. 

When  the  inner  sphere  activation  barrier  is  large,  and/or  when  the  reaction  adiabaticity 
is  weak,  the  dependence  of  the  pre-exponential  factor  on  solvent  is  weakened.  A  convenient 
approximate  way  of  writing  the  nuclear  barrier  crossing  frequency,  ¥n,  is  then  equation  (14) 

[33] 

vn  =  Vis^Ti/8  (14) 

where  is  the  frequency  associated  with  an  intra-molecular  vibration.  Combining 
equations  (8)  and  (14),  one  may  express  the  standard  rate  constant  using  equation  (IS) 
ks  =  AtL-0exp(-AG*/RT)  (15) 

where  the  factor  A  contains  all  the  contributions  to  the  pre-exponential  factor  which  to  a  first 
approximation  do  not  depend  on  the  solvent  and  ft  is  a  fraction  between  0  and  1.  This  way 
of  expressing  the  solvent  dependence  was  first  given  by  McGuire  and  McLendon  [34]  and  is 
especially  convenient  for  discussing  the  experimentally  observed  solvent  dependence.  It  is 
also  suggested  by  the  theoretical  analysis  presented  by  Nadler  and  Marcus  [35].  An 
unattractive  feature  of  equation  (15)  is  that  the  factor  A  has  fractional  dimensions  of  time 
when  ft  is  a  fraction  [33]. 

The  nuclear  barrier  crossing  frequency  on  the  basis  of  known  values  of  Yis  and  II  is 
the  order  of  1012  to  1014  s*1.  As  a  result,  the  upper  limit  for  the  pre-exponential  factor  for 
heterogeneous  ET  is  approximately  6000  cm  s*1.  For  a  reaction  with  a  Gibbs  energy  barrier 
of  20  kJ  mol*1,  this  corresponds  to  a  standard  rate  constant  of  2  cm  s*1  at  25  °C.  This 
estimate  is  certainly  of  the  correct  order  of  magnitude  in  terms  of  kinetic  data  collected  to 
date. 
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Some  mention  should  be  made  of  solvents  which  do  not  follow  the  simple  Debye 
model.  Most  solvents  exhibit  multiple  relaxation  when  the  frequency  dependence  of  the 
dielectric  permittivity  is  examined  over  a  sufficiently  wide  range  [36],  Estimation  of  II  is 
more  complex  for  these  systems  [37],  but  can  be  made  on  the  basis  of  the  work  of  Hynes 
[28].  Solvents  for  which  non-Debye  behaviour  is  important  are  the  alcohols,  and  protic 
amides  [36].  In  these  systems  Jl  is  a  function  of  time  [37],  the  value  observed  at  very  short 
times  (high  frequency  limit)  being  relevant  to  the  type  of  reaction  being  considered  here. 

Thus,  when  II  for  the  alcohols  was  estimated  considering  die  relaxation  modes  at  higher 
frequencies,  the  kinetic  data  obtained  in  these  solvents  falls  into  line  with  those  obtained  in 
simpler  solvents  [33]. 

As  a  note  of  caution,  it  should  be  remembered  that  the  Debye  model  on  which  the 
estimation  of  II  is  based  is  a  continuum  model,  and  therefore  does  not  recognize  the  unique 
properties  of  the  solvent  in  the  vicinity  of  the  reactant  Thus,  the  local  value  of  ^  may  be 
quite  different  from  that  estimated  on  the  basis  of  the  Debye  parameters  (equation  13). 
However,  the  fact  that  one  obtains  clear  correlations  between  In  fcs  and  ]q  for  the  majority  of 

systems  studied  to  date  suggests  that  there  is  a  simple  relation  between  the  actual  local 
relaxation  time  and  that  estimated  by  the  Debye  model.  This  problem  can  be  elucidated  using 
ultra  fast  spectroscopic  relaxation  techniques  and  work  in  that  direction  has  been  initiated  [10]. 

Finally,  it  should  be  noted  that  the  above  discussion  has  been  made  within  the  context 
of  Marcus  theory  which  strictly  speaking  was  developed  for  adiabatic  reactions  [12]. 
Electronic  effects  are  included  in  this  model  heuristically  by  way  of  die  electronic 
transmission  coefficient  &  An  alternative  approach  is  that  of  Dogonadze  [38]  which  is  based 
on  a  quantum  mechanical  description  of  a  manifold  of  potential  energy  surfaces.  The 
quantum  mechanical  approach  was  recently  extended  by  Schmickler  [39]  and  its  relationship 
to  relevant  experimental  observations  discussed  in  detail. 
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3 .  New  Developments  In  Experimental  Methods 

In  order  to  test  the  predictions  made  by  contemporary  theory,  one  must  be  able  to 
determine  the  standard  rate  constants  for  fast  heterogeneous  ET  reactions.  Efforts  in  this 
direction  have  most  often  involved  ultramicroelectrodes  [40-42].  When  these  are  used  as 
working  electrode,  the  current  flowing  in  the  electrochemical  cell  is  typically  100  to  1000 
times  smaller  than  that  at  a  conventional  dropping  mercury  electrode  so  that  problems  due  to 
iR  drop  and  capacitadve  current  are  greatly  reduced.  Ultramicroelectxodes  are  often  metal 
wires  with  diameters  of  the  order  of  microns  which  are  sealed  in  glass  so  that  a  disc  is 
exposed  to  the  electrolyte  solution.  Electrode  materials  include  platinum,  gold  and  carbon, 
and  must  be  carefully  polished  before  each  experiment  The  gold  disc  electrode  may  be 
modified  to  a  mercury  hemisphere  by  electrodeposition  of  mercury  [43,44],  or  by  exposing 
the  gold  disc  to  a  pool  of  highly  purified  mercury  for  a  short  time  [43].  Attempts  have  also 
been  made  to  fabricate  electrodes  of  nanometer  dimensions,  so  called  "nanodes"  [46]. 
However,  the  area  of  the  conducting  material  exposed  to  the  solution  is  extremely  difficult  to 
define  due  to  problems  in  sealing  glass  to  metal  rods  of  such  small  diameter  [47].  As  a  result 
interpretation  of  current-voltage  data  with  a  view  to  obtaining  kinetic  information  is 
impossible  [47,48]. 

The  most  popular  technique  for  determining  kinetic  parameters  using 
ultramicroelectrodes  has  been  cyclic  voltammetry  using  fast  scan  rates  [42, 49, 50]. 

However,  extraction  of  the  kinetic  information  is  not  straightforward  and  requires  careful 
analysis  of  the  effects  of  any  uncompensated  iR  drop  [51]. 

Another  technique  introduced  recently  by  Baranski  [43, 44, 52-54]  is  high  frequency 
a.c.  admittance  voltammetry  at  ultramicroelectrodes.  By  working  at  higher  frequencies  linear 
mass  transfer  conditions  are  maintained  at  the  working  electrode.  Moreover,  die  resistance  of 
the  system,  and  interfacial  capacitance  are  monitored  directly  in  the  experiment  so  that  any 
problems  with  the  state  of  the  electrode/solution  interface  are  easily  recognized.  Baranski 
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[54]  has  carried  out  a  detailed  analysis  of  the  possible  errors  in  these  experiments  and  shown 
how  one  may  optimize  the  frequency  range  for  a  kinetic  study  given  the  electrode  area, 
solution  resistance,  and  double  layer  capacity.  These  experiments  involve  a  fairly  simple 
extension  of  well  known  admittance  methods  [55, 56]  and  are  easily  performed  at 
frequencies  up  to  100  kHz.  As  a  result,  there  is  no  problem  in  determining  rare  constants  up 
to  10  cm  S'1  with  reasonable  precision  [44, 45,  53, 57]. 

Bard  and  coworkers  [58-61]  have  demonstrated  that  scanning  electrochemical 
microscopy  may  be  used  to  determine  fast  ET  kinetics  precisely.  In  such  an  experiment,  a 
micron-sized  ultramicroelectrode,  which  is  normally  used  to  probe  the  surface  of  a  substrate 
and  is  very  close  to  it,  acts  as  a  working  electrode  with  the  reactant  dissolved  in  die  small 
volume  of  solution  between  the  ultramicroelectrode  and  the  closely  positioned  substrate. 

They  have  shown  that  mass  transfer  conditions  are  such  that  one  can  measure  fast  ET  kinetics 
in  both  aqueous  [58, 59]  and  highly  resistive  non-aqueous  media  [60, 61].  Results  from 
these  experiments  for  the  electrooxidation  of  ferrocene  in  acetonitrile  [59]  and  the 
electroreduction  of  buckminsterfullerene  in  benzonitrile  [61]  agree  well  with  data  obtained 
earlier  using  high  frequency  a.c.  admittance  methods  [53, 57]. 

Osteryoung  et  al.  [62, 63]  have  described  a  short  time  pulse  voltammetric  technique 
for  studying  fast  heterogeneous  reactions  at  ultramicroelectrodes.  In  the  normal  pulse  mode, 
the  electrode  potential  is  stepped  to  a  potential  where  a  faradaic  process  occurs  and  held  at  the 
final  potential  for  a  fixed  time  which  is  the  order  of  microseconds,  and  the  resulting  current 
transient  recorded  [62].  The  kinetic  information  is  obtained  from  the  faradaic  current 
measured  towards  the  end  of  the  pulse.  The  electrode  is  held  at  the  initial  potential  for  a  time 
which  is  the  order  of  milliseconds  so  that  the  effects  of  mass  transfer  near  the  electrode  are 
obliterated  by  diffusion  before  the  next  pulse  to  a  new  final  potential  is  applied.  This 
technique  was  recently  applied  to  determining  standard  rate  constants  for  ferrocene  and 
anthracene  in  acetonitrile  and  dichloromethane  [63]. 
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Mahlendorf  and  Heinze  [64]  have  recently  used  impedance  spectroscopy  to  study  fast 
electron  transfer  reactions  at  a  Pt  electrode  of  conventional  size  in  highly  resistive  non- 
aqueous  media.  Their  method  of  data  analysis  differs  significandy  from  that  used  by 
Baranski  [44, 54]  and  requires  the  extraction  of  five  parameters  using  non-linear  regression. 
They  obtained  a  standard  rate  constant  for  the  oxidation  of  ferrocene  in  die hlorome thane 
which  is  slighdy  less  than  that  found  by  a.c.  admittance  techniques  at  ultramicroelectrodes 
[53]  when  one  considers  the  fact  that  their  experiments  were  carried  out  at  a  higher 
temperature  (293  K).  These  authors  [64]  were  unsuccessful  in  obtaining  kinetic  data  in 
acetonitrile. 

Some  feeling  for  the  experimental  problems  involved  in  determining  rale  constants  for 
fast  heterogeneous  reactions  can  be  obtained  by  comparing  data  reported  recently  by  various 
laboratories  for  the  kinetics  of  electrooxidation  of  ferrocene  in  acetonitrile  (see  Table  2). 
Acetonitrile  is  a  solvent  of  low  viscosity  with  a  short  longitudinal  relaxation  time  (~0-2  ps), 
and  ferrocene  is  a  reactant  with  a  low  inner  sphere  reorganization  energy.  These  factors 
combine  to  make  ET  rapid  in  this  system.  Examination  of  die  data  summarized  in  Table  2 
reveals  an  unacceptably  large  range  for  rate  constants  reported  during  the  last  seven  years, 
namely,  from  2xl0'2  cm  s*1  to  220  s*1,  that  is  by  four  orders  of  magnitude!  If  the  two  low 
results  [66, 67]  and  the  high  result  [46]  are  removed  a  much  more  acceptable  range,  from  1 
to  4  cm  s*1,  is  found.  It  should  be  noted  that  differences  in  electrolyte  nature  and 
concentration,  and  in  temperature  cause  some  of  the  observed  scatter  in  the  remaining  data. 
The  very  high  result  by  Penner  et  al.  [46]  is  almost  certainly  due  to  problems  in  defining  die 
area  of  the  electrode  exposed  to  the  solution  when  electrodes  of  diameters  significandy  less 
than  1  pm  are  used  [47, 48, 60],  The  same  comment  can  be  applied  to  the  low  result  of 
Russell  et  al.  [66]  who  used  a  ring  microelectrode.  The  low  result  reported  by  Abbott  et  al. 
[67]  at  an  electrode  of  conventional  dimensions  is  attributed  to  uncompensated  iR  drop.  The 
results  obtained  by  fast  scan  cyclic  voltammetry  at  ultramicroelectrodes  [49, 50, 65]  fall 
somewhat  below  the  most  recent  result  by  Miririn  et  al  [60],  probably  as  a  result  of  some 
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uncompensated  iR  drop  [51]  and  effects  due  to  non-linear  diffusion  [60, 68].  The  rate 
constant  estimated  from  the  data  of  Baranski  et  al.  [53]  on  the  basis  of  available  temperature 
coefficient  data  [53, 69, 70]  is  3.8  cm  s'1  at  298  K.  Thus,  the  standard  rate  constant  for  this 
reaction  determined  by  scanning  electrochemical  microscopy  and  high  frequency  a.c. 
admittance  agree  with  one  another  within  experimental  error.  The  recent  result  of  Karpin  ski 
and  Osteryoung  [63]  falls  somewhat  below  these  for  reasons  which  are  difficult  to  assess 
but  which  also  may  be  related  to  uncompensated  iR  drop. 

The  above  results  demonstrate  that  one  must  exercise  great  care  in  the  measurement  of 
fast  ET  kinetics.  However,  it  is  clear  that  experimental  techniques  are  now  available  which 
can  provide  precise  data  for  rapid  reactions.  This  means  that  important  theoretical  questions 
can  be  examined  on  the  basis  of  reliable  experimental  data 

4.  Effects  Of  The  Solvent 

The  major  focus  of  experimental  work  during  the  last  ten  years  has  bear  to 
understand  the  role  of  the  solvent  [71, 72].  These  studies  have  shown  clearly  that  a  dynamic 
solvent  effect  is  present  especially  when  die  inner  sphere  reorganizations!  energy  is  small. 

The  presence  of  a  dynamic  effect  is  most  easily  demonstrated  by  plotting  die  standard  rate 
constant  measured  in  various  solvents  against  the  solvent’s  longitudinal  relaxation  time  on 
logarithmic  scales.  Data  for  the  electrooxidation  of  ferrocene  [53]  and  1 ,4-diaminobenzene 
[73]  are  shown  in  Figure  1.  Quite  good  linear  relationships  are  obtained  between  these  two 

quantities,  the  slope  for  the  1 ,4-diaminobenzene  system  being  less  than  that  for  the  ferrocene 

* 

system.  The  scatter  observed  is  attributed  to  the  fact  that  die  variation  in  AG^  with  solvent 

is  ignored  in  making  such  a  plot  On  the  basis  of  equation  (15),  the  exact  relationship 
between  In  kc  and  IfllL  is 

In  ks  =  In  A  -  0  In  xL  -  AG*/RT  -  AG^/RT 


(16) 
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The  results  obtained  for  many  systems  [5, 45, 53, 69, 73-79]  support  the  conclusion  that 

* 

variation  in  In  a  is  significantly  greater  than  variation  in  AGQS  /El  fen1  a  typical  selection  of 
polar  solvents.  In  early  work  [5],  it  was  also  pointed  out  that  a  correlation  exists  with  the 
logarithm  of  the  solvent's  viscosity. 

Considerable  discussion  has  centered  on  the  separation  of  die  solvent  effect  due  to  the 

pre-exponential  factor  from  that  due  to  the  Gibbs  energy  barrier  [1 1, 31, 33, 40].  If  both  II 
* 

and  AG0&  are  known,  the  parameters  A  and  ft  are  easily  determined  using  linear  regression 

* 

techniques.  In  most  cases  AG,_  is  either  unknown  or  difficult  to  estimate.  However,  on  the 

basis  of  the  Marcus  model  [12],  it  is  proportional  to  the  solvent's  Pekar  factor,  1/fop  -  !/£«• 

* 

An  approximate  linear  relationship  between  AGQS  and  the  Pekar  factor  is  also  found  when 
the  former  quantity  is  estimated  using  the  FB  model  [22].  Thus,  equation  (16)  can  be  written 
in  the  approximate  form 

lnks  s  K  -01ntL -gf—  *  (17) 

Uop  Es  J 

* 

where  &  is  a  solvent  independent  constant  and  g,  the  proportionality  constant  relating  AGQ£ 

to  the  Pekar  factor.  Thus,  the  parameters  &  g  and  £  can  be  determined  on  the  basis  of  a  two 
parameter  regression  with  la  II  and  the  Pekar  factor  as  independent  variables  [30].  Such  an 
analysis  requires  that  ]qIl  and  Q/Sop  -  i/£s) not  be  correlated  to  one  another,  and  that  data  be 
available  in  a  sufficient  number  of  solvents  (at  least  four).  This  analysis  has  been  applied 
successfully  to  kinetic  data  for  several  heterogeneous  ET  reactions  [30, 45, 77]  and  is 
especially  useful  in  obtaining  an  estimate  of  the  parameter  ft  Weaver  [1 1]  has  strongly 
criticized  this  method  of  analysis,  pointing  out  that  there  is  a  weak  correlation  between  ]q  II 
and  (1/fop  *  1/fo).  This  point  is  illustrated  in  Figure  2  on  the  basis  of  parameters  for  15  polar 
solvents.  A  very  weak  correlation  is  present  for  the  group  of  solvents  considered  (r  « 

0.733),  the  value  of  the  Pekar  factor  becoming  smaller  as  the  value  of]niL  increases. 
Obviously,  any  properly  performed  statistical  analysis  must  take  this  into  consideration,  but 
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the  degree  of  correlation  seen  for  the  IS  solvents  chosen  would  not  preclude  an  analysis  of 
kinetic  data  obtained  provided  the  level  of  experimental  error  in  the  standard  heterogeneous 
ET  rate  constant  is  low.  The  data  presented  in  Figure  2  also  demonstrate  that  one  must  design 
experimental  studies  of  solvent  effects  with  careful  consideration  of  the  degree  of  correlation 
between  ]q  Zl  and  the  Pekar  factor.  For  example,  Phelps  et  al.  [80]  recently  reported  a  study 
of  solvent  effects  on  the  electrooxidation  of  several  sesquibicyclic  hydrazines  in  seven  polar 
solvents,  namely,  AC,  AN,  NM,  DMF,  DMSO,  NB  and  MeOH.  The  correlation  between  la 
II  and  the  Pekar  factor  is  more  important  for  this  group  of  solvents  (r  =  0.89S)  so  that  one 
has  no  hope  of  separating  the  contributions  from  the  pre-exponential  factor  and  the  activation 
barrier  on  the  basis  of  equation  (17)  and  less  objective  methods  must  be  used  to  assess  these 
data.  A  more  careful  selection  of  solvents  used  would  have  allowed  firmer  conclusions  to  be 
made  about  the  role  of  the  solvent  for  these  systems. 

Other  techniques  for  analyzing  the  solvent  effect  involve  assuming  that  one  of  the 

solvent  dependent  terms  is  known,  correcting  the  observed  rate  constant  for  it,  and  then 

examining  the  other  solvent  effect  with  respect  to  the  corrected  rate  constant  [1 1].  For 

* 

instance,  if  one  assumes  values  for  AG  .  one  may  write  equation  (16)  as 

In  ks  +  AGqj/RT  =  In  A  -  AG*  /  RT  -  0  In  tL  (18) 

It  follows  that  a  plot  of  the  left  hand  side  of  equation  (18)  against  la  II  should  yield  the 

parameter  0.  This  procedure  was  followed  in  two  recent  studies  [80, 81]  using  the 

* 

continuum  estimates  of  AGps  (equation  2)  with  £  (imaging)  set  equal  to  infinity.  These 
* 

estimates  of  AGQ£  are  almost  certainly  high  by  at  least  1  £1  and  perhaps  by  as  much  as4RT 

depending  on  the  nature  of  the  reaction,  so  that  the  analyses  are  of  doubtful  value.  A  more 

* 

interesting  method  of  estimating  AGQS  is  by  using  the  optical  energies  for  intramolecular 

electron  transfer  in  bis(ferrocenyl)  acetylene  cation  and  related  compounds  [1 1, 82, 83]. 
Thus,  the  energy  associated  with  transfer  of  one  electron  from  one  ferrocene  group  to  die 
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other  across  the  acetylene  bridge  is  assumed  to  be  directly  related  to  the  value  of  AG*&  for 

homogeneous  ET  between  ferrocene  and  ferrocinium  cation  in  a  given  solvent  The  estimates 
of  AG^  from  these  data  are  approximately  85  percent  of  those  obtained  by  the  continuum 

model  and  therefore  closer  to  estimates  based  on  the  MSA  results  (Table  1).  They  have 

mainly  been  used  in  the  analysis  of  data  for  homogeneous  ET  involving  the  metallocenes 
[11].  This  method  of  estimating  AGQS  is  considered  in  more  detail  below. 

A  method  applied  in  the  case  of  perfectly  adiabatic  reactions  (ft  =  1)  made  use  of 

equations  (2)  and  (12),  so  that  the  standard  rate  constant  was  corrected  for  solvent 

dependence  of  the  pre-exponential  factor  [84, 85].*  *  This  method  is  not  of  great  interest 

now  considering  the  body  of  evidence  which  indicates  that  equation  (2)  overestimates  the 
« 

value  of  AGqs. 

Some  mention  should  be  made  of  the  model  for  solvent  effects  in  heterogeneous 
rate  constants  proposed  by  Abbott  et  al.  [67]  on  the  basis  of  data  for  the  oxidation  of 
ferrocene  in  16  non-aqueous  solvents.  These  authors  ignored  possible  dependence  of  the 
pre-exponential  factor  of  the  rate  constant  mi  solvent  relaxation,  and  attributed  solvent 
dependence  to  variation  in  AG*  with  solvent  basicity,  acidity,  polarity,  and  polarizability. 
Indeed,  the  thermodynamic  properties  related  to  ion  solvation  can  be  related  to  solvent 
basicity  in  the  case  of  cations  and  solvent  acidity  in  the  case  of  anions  [134],  but  ion 
solvation  at  optical  frequencies  has  not  been  considered  on  the  basis  of  the  models  which 
they  discuss.  Considering  the  poor  quality  of  most  correlations  reported  by  Abbott  et  aL  [67] 
and  the  probability  that  their  kinetic  data  are  strongly  influenced  by  iR  drop,  the  conclusions 
reached  in  their  study  are  difficult  to  accept 


*  In  reference  [85],  the  method  was  incorrectly  applied  to  some  reactions  which  are  not  perfectly 
adiabatic  (6  <  1). 


19 


The  general  picture  which  has  emerged  recently  is  that  the  solvent  contribution  to  the 
pre -exponential  factor  can  vary  over  a  wide  range  of  values  of  Q.  Available  data  in  the 
literature  were  examined  by  the  present  authors  [33]  who  showed  that  the  magnitude  of  Q 
was  connected  to  the  magnitude  of  the  rate  constant  observed  in  a  given  solvent 
Accordingly,  as  the  reaction  becomes  slower  the  value  of  Q  decreases  (see  Figure  1).  This  is 
attributed  to  a  larger  contribution  to  the  activation  barrier  from  internal  vibrational  modes, 
and/or  weak  reaction  adiabaticity.  When  $  is  sufficiently  small  one  can  observe  a  linear 
correlation  between  the  logarithm  of  the  standard  rate  constant  and  the  Pekar  factor,  as  noted 
by  Weaver  et  al.  [86]  for  the  tri(hexafluoroacetylacetonato)  ruthenium(III)  system. 

In  conclusion,  although  many  studies  of  solvent  effects  have  now  been  carried  out, 
and  the  presence  of  a  dynamic  effect  clearly  established,  one  cannot  assess  the  role  of  the 
solvent  on  the  basis  of  these  data  alone.  This  is  chiefly  because  one  does  not  have  a  good 
method  of  estimating  the  activation  barrier  except  for  the  simplest  systems.  One  way  around 
this  problem  is  to  determine  the  activation  enthalpy  for  the  reaction  by  studying  the 
temperature  dependence  of  the  heterogeneous  rate  constant  These  experiments  are 
considered  in  the  following  section. 

5.  The  Experimental  Activation  Enthalpy 

Studies  of  the  temperature  dependence  of  the  standard  rate  constant  for  heterogeneous 
ET  are  particularly  useful  when  they  are  available  in  different  polar  solvents,  and  also  for 
both  the  homogeneous  and  heterogeneous  reaction.  On  the  basis  of  equation  (15),  the 
experimental  activation  enthalpy  is  given  by  equation  (19) 

4H«  -  -Rti?T)  “  +  AHos  +  e  Alt  <»> 

* 

where  AH|S  is  the  inner  sphere  activation  enthalpy,  and  AHi .  the  enthalpy  arising  from  die 

temperature  dependence  of  the  longitudinal  relaxation  time.  The  latter  quantity  depends  on 
the  temperature  dependence  of  the  Debye  relaxation  time.  Id.  and  the  high  frequency  and 
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static  dielectric  constants,  £»  and  (equation  13).  Values  of  the  temperature  coefficients  of 
Id  and  g*  are  available  for  most  polar  solvents  [87]  and  that  for  £»  can  be  estimated  from  the 
temperature  coefficient  of  die  refractive  index  of  the  solvent  [88].  AHi  is  not  insignificant 

and  varies  from  ~4  kJ  mol'1  in  fast  solvents  to  -13  kJ  mol*1  in  slow  solvents. 

* 

Solvent  dependent  AHCX  data  have  been  collected  for  a  number  of  heterogeneous  ET 

reactions  including  the  electroreduction  of  cobaltacinium  cation  [69],  cobaltocene  [88],  the 

acetylacetonato  complexes  of  Mnm  ami  Fe111  [77],  benzophenone  [45]  and  a  manganese111 

porphyrin  [79].  Data  are  also  available  for  the  electrooxidation  of  ferrocene  [53, 70]  and 

* 

tetrathiofulvalene  [78].  In  general,  one  finds  that  the  value  of  AHCX  increases  with  increase 

in  solvent  viscosity.  This  result  demonstrates  that  the  dynamic  solvent  effect  expressed 
through  the  AHi  term  in  equation  ( 19)  is  present  in  these  data.  The  most  interesting  system 
from  the  point  of  view  of  heterogeneous  ET  is  cobaltacene  because  this  molecule  may  be  both 
oxidized  to  a  cation  and  reduced  to  an  anion.  On  the  basis  of  the  FB  model,  the  activation 
parameters  for  these  reactions  should  not  be  the  same,  whereas  the  Marcus  model  predicts 
that  they  are  the  same.  In  addition,  activation  enthalpy  data  are  available  for  the 
homogeneous  ET  between  cobaltacinium  cation  and  the  cobaltocene  molecule  [89].  For  these 
reasons  this  system  is  chosen  here  for  more  detailed  examination. 

The  available  data  for  the  redox  reactions  of  cobaltocene  are  summarized  in  Table  3. 

* 

Where  comparisons  can  be  made,  it  is  seen  that  AHCX  for  the  homogeneous  oxidation 

reaction  is  less  than  that  for  the  corresponding  heterogeneous  reaction.  Furthermore,  the 
* 

value  of  AHCX  for  the  heterogeneous  reduction  reaction  is  larger  than  that  for  the  oxidation 

reaction  [88].  The  experimental  results  are  often  less  than  the  continuum  estimates  of*C 

(see  Table  1).  Keeping  in  mind  that  there  are  contributions  to  AH^  from  the  inner  sphere 

activation  enthalpy  and  the  longitudinal  relaxation  enthalpy,  these  data  strongly  confirm  that 

* 

the  continuum  estimates  of  are  too  large.  One  method  of  testing  the  models  considered 

*  *  * 
here  is  to  plot  AHCX  -  AHqS  against  AHi  and  obtain  9  from  the  slope  and  AH^  from  the 
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* 

intercept  [88].  These  plots  using  the  FB  estimates  of  AHQ£  are  shown  in  Figure  3  for  both 
the  homogeneous  and  heterogeneous  data.  Reasonable  linear  plots  are  obtained  for  both 
cases,  die  slope  being  larger  for  the  homogeneous  data.  It  would  seem  that  the  FB  model 
overestimates  the  activation  parameters  for  the  homogeneous  case  perhaps  due  to  effects 

related  to  the  fact  that  the  two  reactants  are  close  to  one  another  in  the  transition  state.  This 

* 

conclusion  is  supported  by  the  observation  that  the  intercept  which  gives  is 
significantly  smaller  than  independent  estimates  of  this  quantity  (2.1  kJ  mol1  [89]).  When 
the  same  analysis  is  applied  to  the  heterogeneous  data,  very  striking  results  are  obtained 
using  the  FB  model,  the  data  for  the  electroreduction  of  cobaltacene  lying  on  die  same  line  as 
those  for  the  electrooxidation  of  the  same  molecule.  These  results  support  the  proposition 
that  differences  in  the  solvation  of  cations  and  anions  at  static  frequencies  must  be  considered 
in  estimating  the  activation  parameters  fen1  ET  reactions  [23]. 


Comparison  of  homogeneous  and  heterogeneous  data  for  AHCX  for  other  systems 

demonstrates  other  features  of  ET  reactions.  In  the  case  of  1 ,4-diaminobenzene,  the  values 
* 

ofAHex  determined  for  the  heterogeneous  reaction  in  AN  and  DMF  [90]  are  much  larger 
than  the  cor responding  quantities  for  the  homogeneous  reaction  [91].  This  is  certainly  due  in 
pan  to  the  fact  that  there  is  a  dynamic  solvent  effect  in  the  heterogeneous  reaction  [73] 
whereas  it  is  absent  for  the  homogeneous  system  [30, 92].  A  more  detailed  analysis  is  not 
possible  because  data  are  available  for  only  two  solvents.  Similar  comments  apply  to  the 
tetrathiafulvalene  system  [78]. 


Particularly  interesting  results  were  described  recently  by  Mu  and  Schultz  [79]  who 

demonstrated  both  on  the  basis  of  rate  constants  and  activation  enthalpies  for  the  reduction  of 

a  Mnm  porphyrin  in  six  aprotic  solvents  that  a  dynamic  solvent  effect  is  present  for  a  system 

*  *  , 

with  a  quite  large  inner  reorganization  energy  (AG-£  =  AHjs  =  16  kJ  tnol*1).  These  results 

suggest  that  factors  other  than  those  usually  considered  [1-4, 9, 1 1, 28]  are  important  in 
determining  the  degree  of  reaction  adiabaticity  and  hence,  the  magnitude  of  the  dynamic 
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solvent  effect.  The  reactant  in  this  study  is  clearly  non-spherical,  and  the  authors  obtained  an 

estimate  of  an  effective  radius  on  the  basis  of  diffusion  coefficient  data.  However,  their 
*  * 

estimates  of  AG0S  and  AH^&  were  based  on  the  continuum  model  and  therefore  are  probably 
too  large.  As  a  result,  the  estimate  of  £  based  on  the  solvent  effect  at  25  °C  was  different 
from  that  obtained  from  the  activation  enthalpy  data  [791.  When  the  FB  model  is  used  to 
estimate  these  quantities,  the  value  of  £  from  the  two  analyses  is  the  same  (-0.7). 

It  is  clear  that  studies  of  the  temperature  dependence  of  the  standard  rate  constant  in  a 
series  of  polar  solvents  gives  an  additional  powerful  tod  for  examining  the  nature  erf  the  ET 
step.  Under  these  circumstances  one  can  assess  the  degree  of  reaction  adiabaticity  provided 
one  has  a  way  of  estimating  the  effective  radius  of  the  reactant  This  subject  is  discussed  in 
the  following  section. 

6 .  Effects  Of  The  Reactant 

When  one  examines  the  data  available  in  the  literature,  one  finds  that  die  standard  rate 

constant  in  a  given  solvent  varies  by  four  orders  of  magnitude  with  the  nature  of  the  reactant 

* 

[33].  Much  of  this  change  is  usually  associated  with  variation  in  AG^  within  the  group  of 

* 

reactants  considered,  but  some  is  also  due  to  changes  in  AG^.  In  order  to  estimate  the  latter 

contribution,  an  estimate  of  the  effective  radius  of  the  reactant  is  needed.  For  this  reason, 

much  experimental  effort  has  been  directed  in  the  recent  past  to  studying  die  metallocenes 

which  are  simple  from  the  point  of  view  that  they  can  be  represented  to  a  good  approximation 

as  spheres  [26].  The  same  comment  applies  to  buckminsterfullerene  whose  heterogeneous 

ET  kinetics  have  been  studied  more  recently  [57, 61, 93).  For  this  system,  quite  good 
*  * 

estimates  of  AG0„  and  AHQ£  were  made  using  the  FB  model  assuming  that  die  effective 
reactant  radius  is  equal  to  the  crystallographic  radius  (360  pm)  plus  the  van  der  Waals  radius 
of  carbon  (70  pm)  [93]. 
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However,  most  reactants  are  not  spherical  and  other  methods  must  be  used  to  deal 
with  the  shape  of  the  molecule  and  the  charge  distribution  within  it.  Early  work  by  Peover 
and  Powell  [94]  involved  the  study  of  the  kinetics  of  electroreduction  of  nitrocompounds. 
As  the  charge  on  the  anion  radical  formed  became  more  delocalized,  say  on  a  benzene  ring, 
the  rate  of  reaction  increased.  A  similar  effect  was  observed  in  the  electroreduction  of  a 
series  of  quinones,  the  standard  rate  constant  increasing  with  increase  in  the  number  of 
benzene  rings  in  the  series  benzoquinone  <  naphthoquinone<anthraquinone  [95]. 


Fawcett  and  Kharkats  [96]  extended  the  model  proposed  by  Peover  and  Powell  [94] 
to  take  into  account  the  volume  correction  terms  when  the  reactant  is  represented  as  two 
contiguous  conducting  spheres.  The  expression  for  the  effective  radius  fc  for  such  a  system 
neglecting  the  effect  of  images  is 

+  fb  +  2  ^  '  f>  F(C-  b>  '  f2  F<C'  »>]  (20> 


where 


F(r*,n) 


(21) 


Here  fi  and  ^  are  the  fractions  of  the  total  charge  on  spheres  1  and  2,  respectively,  g,  the 

radius  of  sphere  1,  fc,  the  radius  of  sphere  2,  and  £,  the  distance  between  the  centers  of  the 

spheres.  For  contiguous  spheres,  it  is  obvious  that  £  is  equal  to  g  +  fc.  The  first  three  terms 

represent  the  effect  of  the  charges  on  the  spheres  without  consideration  of  the  fact  that  part  of 

the  solvation  atmosphere  of  a  given  sphere  is  excluded  by  the  proximity  of  the  other  sphere. 

The  fourth  and  fifth  terms  account  for  the  effects  of  the  finite  volume  of  each  sphere.  It  is 

clear  that  the  excluded  volume  terms  result  in  an  increase  in  the  effective  radius  of  the 

* 

reactant,  and  thus,  in  a  decrease  in  AG  The  effect  is  small  if  the  radii  of  the  two  spheres 


are  equal  but  increases  in  importance  as  the  relative  sizes  of  die  spheres  changes  [96].  The 
effective  radius  can  be  used  in  the  MSA  model  to  obtain  an  improved  estimate  of  AGQ&.  for 


example,  by  replacing  ia  in  equation  (7)  by  Ce- 
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Fawcett  and  Fedurco  [45]  recently  extended  the  above  model  to  the  benzophenone 
system  which  was  represented  as  three  contiguous  spheres  with  each  sphere  in  contact  with 
two  others.  Two  equal  size  spheres  enclosed  each  of  the  phenyl  rings  and  a  smaller  sphere 
surrounded  the  carbonyl  group.  In  this  case  the  excluded  volume  effect  is  more  important. 
When  the  diameters  of  the  contiguous  spheres  are  defined  by  the  structure  of  the  molecule, 
such  a  model  probably  underestimates  the  extent  of  the  solvation  because  of  the  chosen 
spherical  shape  for  each  part  of  the  molecule.  A  better  representation  would  involve 
contiguous  ellipsoids,  each  extending  somewhat  further  beyond  die  extreme  parts  of  the 
molecule.  However,  a  model  with  ellipsoids  would  involve  more  adjustable  parameters  and 
therefore  would  be  less  easily  justified  by  experiment 


The  excluded  volume  effect  can  also  be  used  to  explain  the  fact  that  the  activation 

enthalpy  for  the  homogeneous  ET  involving  the  cobaltacinium  (+/0)  system  is  slightly 

smaller  than  that  for  heterogeneous  ET  (see  Table  3).  On  the  basis  of  the  woritof  Kharkats 
* 

[97],  AQqj  for  homogeneous  transfer  of  one  electron  between  a  spherical  reactant  of  radius. 


and  another  radius,  fe,  is  given  within  the  context  of  Marcus  continuum  theory  by 
(equation  22) 


AG*  =  NqCq  (—  - 
os  l6jteoTe  tap 

where 

re  =  [l  +  S  -  |-F(c,a)-F(c,b)] 


(22) 


(23) 


The  quantities  £  and  £(cx,  c)  have  the  same  meaning  as  given  above.  When  the  spheres  have 

* 

equal  radii  and  are  in  contact  the  value  of  AGps  is  six  petcent  smaller  than  predicted  when 

the  excluded  volume  terms  are  ignored.  This  correction  can  also  be  made  within  the  context 
of  the  MSA,  the  resulting  relationship  being 
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A  different  approach  to  dealing  with  the  effects  of  molecular  shape  has  been  used  by 

Grampp  et  al.  [78, 91, 92, 98].  In  their  work,  the  molecule  is  represented  as  an  ellipsoid 

with  semi  axes  a,  £  and  £.  This  approach  has  been  applied  especially  to  data  for 

homogeneous  self-exchange  reactions  involving  organic  molecules  such  as 

1,4-diaminobenzene  and  its  derivatives,  tetracyanoethylene,  tetracyano-p-quinonedimethane, 

and  tetrathiofulvalene  [98]  but  also  to  heterogeneous  ET  reactions  [78].  Such  a  description 

of  the  system  requires  more  parameters  than  a  model  based  on  two  continguous  spheres. 

When  applied  to  a  diabatic  system  with  data  available  both  as  a  function  of  solvent  and 

temperature  [91, 92],  a  detailed  assessment  of  the  mutual  orientation  of  the  reactants  in  the 

transition  state  is  possible  [98].  To  date  this  approach  has  been  limited  to  the  Marcus 

* 

continuum  model  for  AGQS  but  it  could  be  easily  extended  to  a  non-primitive  model  such  as 
the  MSA. 


It  should  be  emphasized  that  the  above  models  have  treated  the  reacting  system  as  a 

conductor.  There  may  be  circumstances  when  such  an  assumption  is  not  appropriate  as 

discussed  by  German  and  Kuznetsov  [99, 100].  These  authors  have  considered  the  effects 

* 

of  polarizability  on  the  estimation  of  AGQS ,  and  also  developed  models  in  which  the  reactant 

is  represented  as  a  dielectric  cavity  of  given  shape  with  embedded  point  charges.  Hupp  et  aL 
[101]  have  shown  that  a  dielectric  cavity  model  is  successful  in  accounting  for  the  solvent 
dependence  of  optical  intervalence  charge-transfer  absorption  frequencies  for  symmetrical 
mixed-valence  compounds.  Charge  transfer  within  the  species  (NH3>5Ruffl-4,4,-bipyridine- 
Run(NH3)55+  was  treated  as  dipole  inversion  within  a  low  dielectric  cavity,  the  length  of  the 
dipole  (0.51  nm)  being  based  on  intervalence  Stark  effect  measurements  [102].  On  the  other 
hand,  a  model  based  on  conducting  spheres  separated  by  a  fixed  distance  (equation  22)  did 
not  give  a  reasonable  description  of  the  same  data.  These  authors  obtained  their 
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spectroscopic  data  in  very  dilute  solutions  thereby  avoiding  the  effects  of  ion  pairing  on  the 
spectral  features.  If  the  model  based  on  two  conducting  spheres  is  valid,  the  absorption 
frequency  observed  in  various  solvents  should  be  linear  in  the  solvents'  Pekar  factor  with  a 
slope  determined  by  the  geometrical  parameters  for  the  system  (equation  23).  Although  this 
plot  is  linear  [101]  the  slope  is  much  less  than  what  one  would  expect  on  the  basis  of 
reasonable  values  for  the  size  parameters.  A  similar  result  is  obtained  on  the  basts  of  data  for 
the  bis(ferrocenyl)  acetylene  cation  [11, 82, 83].  The  analysis  carried  out  by  Hupp  et  al. 

[101]  is  quite  important  because  it  suggests  that  model  compounds  based  on  a  redox  couple 
connected  through  a  non-conducting  bridge  are  not  good  systems  for  estimating  the  Gibbs 
activation  energy  for  homogeneous  ET  involving  the  same  redox  couple. 

Other  compounds  for  which  interesting  structural  effects  are  seen  in  association  with 
ET  are  metalloporphyrins  [79, 134-137]  and  oligo-decker  sandwich  compounds  [138].  An 
important  question  in  metalloporphyrin  redox  chemistry  is  whether  ET  takes  place  at  the  metal 
ion  or  to  the  porphyrin  ring.  Mu  and  Schultz  [133]  showed  that  the  porphyrin  centered 
oxidation  of  Mnm  tetraphenylporphyrin  chloride  is  much  faster  than  the  metal  centered 
reduction  of  the  same  compound.  This  observation  was  attributed  to  the  fact  that  the  effective 

radius  of  the  reactant  is  much  larger  in  the  oxidation  reaction  in  which  the  charge  density  is 
distributed  over  the  porphyrin  ring  and  therefore,  AG^S  correspondingly  smaller. 

Anxolab6h6re  et  al.  [136]  obtained  kinetic  data  for  the  reduction  of  some 
tetraphenylporphyrins  and  related  compounds  of  more  complex  structure  in  DMF,  and 
concluded  that  solvent  reorganization  plays  an  important  role  in  these  fast  heterogeneous  ET 
reactions.  Fawcett  et  al.  [  1 37]  examined  the  redox  chemistry  of  a  stable  Zn  isoporphyrin,  and 
concluded  on  the  basis  of  the  kinetic  parameters  and  their  comparison  with  those  for  Zn 
octaethylporphyrin  that  ET  takes  place  on  the  ring  rather  than  to  die  central  metal  km. 
Interpretation  of  these  data  is  greatly  facilitated  when  activation  parameters  are  measured  [79, 
137]. 
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Michels  and  Schmickler  [138]  have  reported  a  study  of  heterogeneous  ET  involving 

oligo-decker  sandwich  compounds.  These  molecules  are  based  on  metallocenes  with  2,3* 

dihydro- 1,3-diborol  as  a  bridging  ligand,  a  basic  example  being  1 ,3,4,5-tetraethyl-2-methyl- 

l,3-diborolyl-bis(cyclopentadienyl)  iron  (BIFc).  Many  of  these  molecules  can  be  both 

oxidized  and  reduced,  and  undergo  some  degree  of  inner  sphere  reorganization  [138].  In  the 

case  of  the  iron  compound  (BIFc),  the  rate  constant  for  the  oxidation  process  was  20  times 

larger  than  that  for  reduction  at  an  Au  electrode  in  die hlorome thane.  The  corresponding 

activation  enthalpies  were  13.3  and  18.3  kJ  mol'1  for  the  oxidation  and  reduction  steps, 

respectively.  These  results  suggest  that  inner  sphere  reorganization  is  not  important  for  the 
iron  compound,  and  support  the  general  conclusion  of  the  FB  model  that  AG^  to  form  a 

cation  in  oxidation  of  a  molecule  should  be  less  than  that  to  form  an  anion.  Further  studies  of 
these  systems  would  be  interesting  with  respect  to  this  question  and  others  discussed  in  this 
review. 


In  summary,  considerable  progress  has  been  made  in  recent  years  to  develop  reliable 
methods  of  estimating  the  effective  size  of  the  reactant  involved  in  heterogeneous  ET.  As  die 
theory  of  ion  solvation  develops  to  deal  with  non-sphcrical  ions  on  the  basis  of  non-primitive 
models  one  can  expect  a  corresponding  improvement  in  understanding  die  role  of  reactant 
shape  in  determining  ET  parameters. 

7.  Double  Layer  Effects 

Most  recent  studies  of  heterogeneous  ET  have  involved  redox  couples  in  which  die 
most  highly  charged  species  is  a  monovalent  cation  or  anion.  Thus,  an  organic  or 
organometallic  molecule  is  involved  in  the  reaction.  This  choice  is  made  in  order  to  minimize 
double  layer  effects  and  effects  due  to  ion  pairing  in  the  bulk  of  the  electrolyte  solution.  For 
this  reason  much  less  attention  has  been  paid  to  redox  couples  involving  transition  metal 
complexes  of  higher  charge  in  spite  of  the  fact  that  such  systems  are  often  easier  to  model 
with  respect  to  structural  effects. 
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The  choice  of  electrolyte  is  also  important  in  establishing  electrostatic  conditions  at  the 
electrode/solution  interface.  In  organic  solutions,  the  electrolyte  is  most  often  a 
tetraalkylammonium  (TAA+)  salt  with  a  hard  anion  such  as  perchlorate,  tetrafluoroborate  or 
hexafluorophosphate.  These  systems  are  chosen  to  avoid  ion  pairing  problems  which  occur 
when  smaller  ions  are  used,  or  problems  with  ionic  adsorption  on  the  electrode.  For 
example,  the  anion  radicals  formed  from  nitro  compounds  strongly  associate  with  alkali  metal 
cations  [103, 104].  In  addition,  the  presence  of  halide  ions  in  the  electrolyte  makes  the 
interpretation  of  kinetic  data  obtained  positive  of  the  potential  of  zero  charge  very  difficult 
because  these  ions  are  strongly  adsorbed  on  polarizable  metals  from  all  known  polar  solvents 
[  105].  TAA+  ions  also  adsorb  on  the  electrode  but  the  nature  of  the  adsorption  appears  to  be 
more  electrostatic  rather  than  chemical  in  character  with  solvophobicity  also  playing  a  role. 
These  ions  also  form  aggregates  with  some  anion  radicals,  notable  examples  being  the  strong 
ion  associates  found  between  C«)'  and  Car  anions  and  TAA+  cations  [57]. 

When  heterogeneous  ET  occurs  at  potentials  negative  of  the  point  of  zero  charge  in 
the  presence  of  TAA+  ions,  it  is  found  that  the  standard  rate  constant  decreases  with  increase 
in  the  size  of  the  alkyl  group  [106-113].  On  the  other  hand,  the  activation  enthalpy  is 
relatively  independent  of  cation  size  [45, 106, 108, 1 12, 113].  Another  observation  is  that 
the  relative  magnitudes  of  the  inhibitory  effect  increases  as  the  electrode  potential  is  made 
more  negative  [1 1 1].  These  observations  can  be  rationalized  if  it  is  assumed  that  the  TAA+ 
ions  are  contact  adsorbed  on  the  electrode  to  form  a  monolayer  whose  thickness  and  extent 
depend  on  the  size  of  the  ion  and  electrode  potential  [111-1 13].  On  the  basis  of  capacity 
data,  it  was  suggested  [113]  that  the  alkyl  arms  of  these  ions  are  bent  when  the  arm  contains 
more  than  three  carbon  atoms  so  that  the  distance  of  closest  approach  of  the  central  nitrogen 
atom  is  approximately  constant  for  the  larger  ions.  Using  this  model,  it  was  possible  to 
explain  the  decrease  in  standard  rate  constant  on  the  basis  of  electron  tunnelling  through  a 
layer  of  TAA+  ions  [114-116].  The  fact  that  the  activation  enthalpy  depends  only  on  the 
solvent  and  not  on  the  size  of  the  TAA+  ion  demonstrates  that  ET  takes  place  outside  of  die 
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adsorbed  ionic  layer  in  a  region  where  die  reactant  and  product  are  still  interacting  chiefly 
with  surrounding  solvent  molecules  [111,113]. 

Other  recent  work  on  electron  tunnelling  through  barrier  films  has  involved  self- 
assembled  monolayers.  Miller  et  al.  [1 17, 1 18]  used  CD-hydroxy  thiols  to  form  a  self- 
assembled  monolayer  on  a  gold  electrode.  The  thiol  functionality  at  one  end  of  he 
hydrocarbon  chain  reacts  with  gold  to  form  a  mercaptide  thereby  anchoring  the  molecule  to 
the  surface.  Interactions  between  the  hydroxy  groups  at  die  other  end  of  die  molecule  and 
water  serve  to  make  the  monolayer  compact  These  authors  [117, 118]  showed  that  electron 
tunnelling  is  the  main  mechanism  for  ET  to  a  redox  system  outside  of  the  monolayer  in 
solution.  By  varying  the  chain  length  of  the  molecule  forming  the  self-assembled  monolayer 
they  were  able  to  demonstrate  that  the  tunnelling  current  depends  exponentially  on  film 
thickness.  This  corresponds  to  an  exponential  decrease  in  the  electronic  transmission 
coefficient  k  in  equation  (8)  [1 15, 1 16]. 

Chidsey  [  1 19]  used  a  self-assembled  monolayer  on  gold  with  an  electroactive 

terminal  group  on  a  fraction  of  the  anchored  groups  to  study  ET  through  films  of  varying 

thickness.  The  electroactive  terminal  group  was  ferrocene,  which  undergoes  fast  ET  upon 

oxidation  to  ferrocinium  cation.  The  effect  of  the  monolayer  consisting  of  hydrocarbon 

chains  with  15  CH2  groups  was  to  slow  the  reaction  down  sufficiently  that  ET  kinetic  data 

could  be  obtained  over  a  voltage  range  of  almost  2  V.  Analysis  of  the  dependence  of  the  rate 

constant  on  potential  provided  evidence  that  the  Gibbs  activation  energy  is  indeed  a  quadratic 

function  of  potential  as  predicted  by  Marcus  theory  [12].  Chidsey  estimated  the  experimental 

value  of  AG*  for  the  ferrocene/ferrocinium  system  in  water  to  be  20.5  kJ  mol'1.  This  result 

* 

is  significantly  lower  than  the  estimate  of  AGQS  in  water  by  continuum  theory  (25.8  kJ  mol* 

1).  The  estimate  of  the  same  quantity  using  the  FB  model  is  17.4  kJ  mol'1.  Considering  the 
fact  that  a  small  contribution  from  inner  sphere  reorganization  must  be  added  (~1  kJ  mol*1), 
the  experimental  result  is  very  close  to  the  value  predicted  by  the  MSA. 
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In  a  similar  study  by  Finklea  and  Hanshew  [120],  the  redox  group  was  pentamine 
(pyridine)  Ru3+/2+.  They  found  some  variation  in  the  experimental  estimate  of  AG* 
depending  on  the  length  of  the  hydrocarbon  chain  and  whether  or  not  spacer  chains  were 

present  in  the  self-assembled  monolayer.  This  suggests  that  the  role  of  the  low  dielectric 

monolayer  in  determining  the  value  of  AG^S  is  changing  in  the  systems  studied  Their 

* 

results  also  support  the  conclusion  that  AGQS  is  overestimated  by  the  continuum  model 
(equation  2). 

Smith  and  White  [121]  have  discussed  double  layer  effects  on  both  the 
thermodynamic  and  kinetic  properties  of  a  self-assembled  monolayer  modified  with  a  redox 
couple.  Since  the  redox  couple  is  not  present  in  the  bulk  of  the  solution,  the  equilibrium 
potential  for  the  surface  attached  redox  couple  depends  on  the  local  (micro)  potentials  at  the 
oxidized  and  reduced  forms  of  the  redox  couple  when  these  are  significantly  different  from 
the  average  potential  in  the  bulk  of  the  solution.  These  authors  [121]  estimated  the  avenge 
potential  on  the  plane  through  the  centers  of  the  redox  species  &  on  the  basis  of  measurable 
thermodynamic  quantities,  and  showed  that  when  this  potential  is  used  to  approximate  the 
local  potential,  a  reasonable  description  of  the  variation  in  equilibrium  potential  of  the  redox 
couple  with  dielectric  properties  of  the  interface  is  possible.  They  also  considered  the  effect 
of  the  charge  associated  with  the  redox  couple  on  the  current  observed  in  linear  sweep 
voltammetry.  The  shape  of  the  current-voltage  curve  was  shown  to  depend  on  the  dielectric 
properties  of  the  interfacial  region,  their  potential  dependence,  and  the  potential  of  zero 
charge.  The  calculations  performed  were  based  on  the  classical  Frumltin  theory  [122, 123] 
for  double  layer  effects  in  electrode  kinetics.  Creager  and  Weber  [124]  extended  this  work 
by  considering  double  layer  effects  on  the  current-time  transients  observed  in  potential  step 
experiments  involving  a  self-assembled  monolayer  modified  with  a  redox  couple. 

The  Frumkin  equation  for  the  logarithm  of  the  double  layer  corrected  forward  rate 
constant  of  the  redox  couple  given  in  equation  (1)  is  usually  written  by  equation  (25) 
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In  kf  +  z*f$d  =  In  ks  +  Ogf^  -  E  -  Ej)  (25) 

where  )tf  is  the  forward  rate  constant,  #*,  the  average  potential  drop  across  the  diffuse  layer, 
2a,  the  charge  on  the  reactant,  E,  the  electrode  potential,  E&>  the  standard  electrode  potential, 
fla,  the  apparent  transfer  coefficient,  and  f  =  F/RT.  The  transfer  coefficient  is  termed 
apparent  because  the  actual  potential  experienced  by  the  reactant  may  differ  from  for  a 
variety  of  reasons  which  arc  generally  termed  discreteness-of-charge  effects  [125-127]. 

These  effects  are  most  noticeable  when  discrete  charges  are  fixed  at  a  dielectric  discontinuity 
in  the  system  [128]  such  as  the  electrode/solution  interface  for  an  unmodified  electrode,  or 
the  hydiocarbon/solution  interface  in  the  systems  studied  by  Chidsey  [1 19]  and  Finklea  and 
Hanshew  [120].  As  a  result  the  local  potential  at  an  anchored  redox  system  jgp may  differ 
significantly  from  the  average  value  on  the  plane  through  their  centers,  $p.  When  this  is  true, 
the  potential,  $d,  in  equation  (25)  should  be  replaced  by  the  local  potential  ]gp  [127]. 

Becka  and  Miller  [  129]  studied  double  layer  effects  at  uncharged  self-assembled 
monolayers  and  showed  that  the  kinetics  of  heterogeneous  ET  at  the  monolayer/solution 
interface  depends  both  on  the  charge  on  the  reactant  and  the  composition  of  the  monolayer. 
The  reduction  of  Fe(CN)g-3  at  a  co-hydroxyl  dodecane  thiol  monolayer  was  considerably 
faster  than  at  a  dodecanethiol  monolayer  at  potentials  negative  of  the  potential  of  zero  charge. 
The  authors  attributed  the  difference  to  the  fact  that  the  self-assembled  monolayer  terminated 
with  the  hydroxy  group  is  less  hydrophobic  so  that  the  reactant  can  approach  closer  to  the 
insulating  film.  In  the  case  of  Fe(bipyridyl)CN4',  the  opposite  effect  is  observed,  die 
reaction  being  faster  at  the  methyl  terminated  monolayer.  This  reactant  is  more  hydrophobic 
and  bears  a  smaller  charge.  These  results  can  also  be  partially  attributed  to  discreteness-of- 
charge  effects.  Because  of  the  large  change  in  dielectric  constant  from  a  low  value  in  die  oil¬ 
like  monolayer  (~2 )  to  a  high  value  in  the  aqueous  solution  (-78),  the  reactant  experiences  a 
strong  repulsive  effect  at  the  monolayer/solution  boundary.  The  electrostatic  effect  [128]  is 
much  more  important  for  the  polyvalent  Fe(CN)63'  ion. 
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Discreceness-of-chargc  effects  have  also  been  discussed  for  heterogeneous  ET 
reactions  occurring  at  mercury  electrodes  in  the  presence  of  specifically  adsorbed  non¬ 
reacting  ions  [130-133].  Three  of  the  systems  studied  [131-133]  involved  transition  metal 
complexes  in  which  the  metal  ion  is  reduced  from  the  +3  to  +2  oxidation  state  in  the  presence 
of  specifically  adsorbed  anions.  By  carrying  out  these  experiments  at  constant  ionic  strength, 
and  examining  the  changes  in  the  logarithm  of  the  corrected  rate  constant,  y  =  inkf + ZAftt*. 
with  adsorbed  charge  density,  fla.  at  constant  potential  drop  across  the  inner  layer,  £  -  £j  - 
$4,  one  can  examine  the  electrostatic  effects  of  the  adsorbed  ions  with  a  minimum  in  change 
in  other  double  layer  properties  (see  equation  25).  The  fact  that  the  adsorbed  ions  are  located 
at  the  metal/solution  interface  means  that  their  electrostatic  accelerating  effect  is  significantly 
reduced  by  the  shielding  effect  of  the  metal.  If  the  simple  Frumkin  equation  were  correct,  the 
kinetic  parameter,  y,  would  be  independent  of  £  -  £$  -  ^  for  changes  in  the  surface  coverage 
by  adsorbed  anions.  In  the  case  of  CofNItyF2*,  y  decreases  with  increase  in  &  as 
qualitatively  expected  on  the  basis  of  the  discreteness-of-charge  effect  [133].  On  the  other 
hand,  when  an  anion  such  as  tetrathionate  is  reduced  in  the  presence  of  adsorbed  iodide 
anions,  the  decelerating  effect  of  the  adsorbed  charge  is  overestimated  when  discreteness-of- 
charge  effects  are  ignored,  and  the  value  of  y  increases  with  increase  in  fl,  [130].  Further 
analysis  of  these  data  [  1 30- 1 33]  led  to  a  more  detailed  description  of  the  double  layer  effect 
including  the  location  of  the  reaction  site  with  respect  to  adsorbed  ions  and  counter  ions. 

Double  layer  effects  can  be  quite  important  in  heterogeneous  ET  reactions  especially 
when  the  reactants  are  highly  charged.  Analysis  of  these  effects  requires  that  interfacial 
thermodynamic  data  be  available  for  the  electrolyte  system  being  used.  The  lack  of  these  data 
have  often  been  the  reason  why  no  attempt  has  been  made  to  examine  the  double  layer  effect 
in  more  detail.  However,  modem  experimental  techniques,  especially  those  based  on 
impedance  allow  one  to  obtain  these  data  relatively  easily.  Double  layer  effects  at  self- 
assembled  monolayers  are  just  beginning  to  be  investigated  and  provide  excellent  systems  on 
which  theories  which  have  been  in  the  literature  for  some  time  can  be  tested  in  more  detail. 
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8 .  Conclusions 

Studies  of  solvent  effects  for  heterogeneous  ET  reactions  have  done  much  to  increase 
the  understanding  of  fundamental  aspects  of  this  elementary  chemical  reaction.  This  is 
especially  true  when  kinetic  data  are  obtained  both  as  a  function  of  solvent  and  of  temperature 
so  that  the  variation  in  activation  enthalpy  may  also  be  examined.  Studies  of  die  redox 
reactions  of  the  metallocenes  [44, 53,  69, 70, 74,  88, 89, 1 10]  have  been  especially 
important  because  of  the  fact  that  the  reactants  are  relatively  spherical  in  shape  and  undergo 
very  small  internal  reorganization  on  ET. 

There  is  strong  evidence  on  the  basis  of  recent  experimental  work  that  the  activation 
parameters  for  outer  sphere  reorganization  are  overestimated  on  the  basis  of  die  original 
Marcus  formulation  [12].  This  is  not  at  all  surprising  since  it  is  well  known  from  equilibrium 
solvation  data  few  simple  ions  that  the  Bom  model  gives  unrealistically  high  estimates  of  the 
Gibbs  energy  and  enthalpy  of  solvation.  The  most  important  feature  of  any  improvement  on 
the  Bom  model  is  that  it  recognizes  the  differences  in  the  solvation  of  cations  and  anions. 
Correspondingly,  one  expects  the  rate  constant  for  formation  of  a  cation  from  a  molecule  to 
be  faster  than  that  for  formation  of  an  anion,  all  other  factors  remaining  the  same.  The 
observation  that  this  is  the  case  for  cobaltacene  [88, 1 10]  and  related  compounds  [138]  is 
important  This  prediction  should  be  examined  for  other  molecules  for  which  these  reactions 
can  be  carried  out. 

Examination  of  the  recent  literature  reveals  that  heterogeneous  ET  can  be  understood 
in  at  least  as  much  detail  as  homogeneous  ET.  Although  one  must  consider  double  layer 
effects  in  the  former  case,  these  can  be  reduced  in  importance  by  careful  choice  of  the 
electrolyte.  Homogeneous  ET  kinetics  has  die  additional  complication  that  the  reactants  may 
assume  a  particular  configuration  with  respect  to  one  another  in  the  transition  state  [98]. 

Such  steric  effects  need  to  be  examined  in  more  detail,  but  they  should  be  simpler  in  die 
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heterogeneous  ET  reaction  where  the  electrical  field  due  to  the  charge  on  the  electrode  plays 
an  important  role  [106].  Changes  in  the  nature  of  the  scene  effect  with  the  type  of  ET 
reaction  may  be  the  cause  of  the  differences  in  the  relative  ordering  of  rate  constants  for 
various  metallocenes  seen  in  the  homogeneous  and  heterogeneous  regimes  [11]. 

Self-assembled  monolayers  with  anchored  redox  groups  [1 19, 120]  has  provided  an 
interesting  new  way  of  studying  heterogeneous  ET.  In  these  systems,  the  environment  of  the 
redox  reaction  is  entirely  different  than  that  at  the  electrode/solution  interface.  Thus,  die 
comparison  of  data  obtained  in  the  two  systems  will  provide  valuable  information  on  the  role 
of  environmental  effects  in  heterogeneous  ET. 

In  conclusion,  the  understanding  of  ET  has  advanced  considerably  during  the  last 
decade,  and  provided  a  fruitful  area  for  both  theoretical  and  experimental  work.  However, 
many  problems  remain  to  be  solved,  so  that  work  in  this  area  will  continue  to  present  a 
challenge  to  those  who  are  involved. 
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Table  1.  Gibbs  energy  and  enthalpy  of  activation  for  a  simple  heterogeneous  electron  transfer  reaction  in  selected  polar 

solvents  according  to  the  Marcus  [12],  McManis- Weaver  (MW)  [20],  and  Fawcett-Blum  (FB)  models  [23).  The 

$  t 

reactant  radius  is  assumed  to  be  0.37  nm  and  the  imaging  distance,  «>.  The  units  for  AGqs  and  AHqs  are  kJ  mol*1. 
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Table  2.  Kinetic  Parameters  for  the  Electrooxidation  of  Ferrocene  at  a  Platinum  Electrode  in  Acetonitrile  from  Recent 
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[a]  TEAP  and  TBAP  are  abbreviations  for  tetraethylammonium  perchlorate  and  tetrabutylammonium  perchlorate,  respectively.  TBABF4 
refers  to  tetrabutylammonium  tetrafluoroborate,  and  TBAPFg  to  tetrabutylammonium  hexafluorophosphate. 


Table  3.  Experimental  Values  of  the  Activation  Enthalpy  AHCX  for  Various  Redox  Reactions  of  Cobaltacene  in  Several 
Aptotic  Solvents  Together  with  the  Longitudinal  Relaxation  Enthalpy  for  the  Solvent  AHl 
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«* 


Figure  Legends 

Figure  I.  Plots  of  the  standard  rate  constant  for  the  electrooxidation  of  ferrocene  (o) 

[53]  and  1 ,4-diaminobenzene  (♦)  [73]  obtained  in  various  aprodc  solvents 
against  the  longitudinal  relaxation  time  of  the  solvent  using  logarithmic  scales 
for  both  axes. 

Figure  2.  Plot  of  the  Pekar  factor  for  a  polar  solvent,  y,  against  its  longitudinal 

relaxation  time  tl  on  a  logarithmic  scale.  The  abbreviations  for  die  solvents 
are  given  in  Table  1  except  for  ethanol  (EtOH),  N-methylpynolidone  (NMP) 
and  proprionitrile  (PN). 

Figure  3.  Plot  of  the  difference  between  the  experimental  activation  enthalpy  AH^  and 

* 

the  estimate  of  the  outer  sphere  contribution  according  to  the  FB  model  AHos 
for  redox  reactions  involving  cobaltacene  against  the  activation  enthalpy  for 
longitudinal  relaxation  AHl-  The  homogeneous  reaction  is  oxidation  of 
cobaltacene  (o)  and  the  heterogeneous  reaction  both  oxidation  (A)  and 
reduction  (V).  The  solid  lines  were  obtained  by  linear  regression  using  all  die 
homogeneous  (o)  or  heterogeneous  data  (A,  V).  The  ordinate  scale  for  the 
heterogeneous  data  has  been  shifted  vertically  by  10  kJ  mol'1  for  clarity.  The 
abbreviations  for  the  solvents  are  given  in  Table  1. 
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